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ABSTRACT 

The  present  work  is  an  extension  of  previous  work  on  the  statistics  of 
quadratic  nonliner  ship  response  to  short-crested  seas.  In  this  previous  work 
the  basic  response  model  was  a  truncated  functional  series;  the  short-crested 
exciting  seaway  was  modeled  is  the  superposition  of  linear  component 
waves.  To  be  mathematically  consistent  the  excitation  must  be  allowed  the 
same  degree  of  nonlinearity  as  the  system  being  excited.  A  consistent 
approach  is  desirable,  even  if  it  merely  yields  indirect  estimates  of  the 
magnitude  of  the  error  resulting  from  the  inconsistency. 

The  objectives  of  the  present  work  included  the  development  of  a 
model  for  the  statistics  of  quadratic  ship  system  response  in  short-crested 
seas  which  is  consistent  in  that  quadratic  nonlinearities  in  the  wave  field  are 
taken  into  account.  It  was  found  that,  if  the  statistical  representation  of  the 
response  to  the  nonlinear  short-crested  deep  water  wave  field  is  carried  out 
strictly  through  contributions  which  may  be  considered  second  order  in 
wave  excitation,  the  statistical  estimating  formulae  for  the  quadratic 
response  to  short-crested  seas  in  which  quadratic  nonlinearities  are  included 
are  identical  to  those  previously  developed  when  no  wave  nonlinearities 
were  included.  The  difference  is  that  the  sum  and  difference  frequency 
response  functions  required  by  the  previous  work  are  modified  by  a  term 
which  amounts  to  the  linear  ship  system  response  to  those  quadratic  wave 
field  components  which  are  produced  by  the  wave  field  nonlinearities.  In 
practice  the  most  important  statistic  involving  quadratic  response  is  the 
expectation,  or  mean  value.  In  this  special  case  an  extension  is  not 
necessary.  To  second  order,  the  mean  response  of  a  linear  plus  quadratic 
ship  system  to  a  short-crested  wave  field  is  not  affected  by  wave  field 
nonlinearities.  The  present  development  suggests  that  the  mean  value 
estimator  may  be  robust  through  third  order. 
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This  work  was  sponsored  by  the  General  Hydromechanics  Research  program  under  Task 
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INTRODUCTION 

Since  the  first  introduction  to  the  marine  field,  [1],  [2]  there  have  been  developed  a 
number  of  applications  of  the  quadratic  functional  model  to  the  prediction  or  simulation  of 
nonlinear  response  of  ships  and  platforms  to  random  seas.  The  details  of  the  various 
approaches  to  quadratic  nonlinear  response  vary.  However,  all  the  state-of-the-art  models  now 
in  use  are  essentially  single  input  models,  meaning  that  validity  is  confined  to  the  long-crested 
case. 

The  basic  formulation  for  the  short-crested  case  was  outlined  in  1966  by  Hasselmann  [1]. 
Recently  this  formulation  was  extended  and  embellished  to  an  extent  by  the  present  writer  [3]. 
This  short-crested  model  involves  the  assumption  that  the  sea  is  a  zero-mean  stationary  and 
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homogeneous  Gaussian  process,  which  excites  a  nonlinear  system  containing  quadratic  level 
nonlinearities  and  thus  produces  a  response  which  is  neither  zero  mean  nor  Gaussian.  As  in  the 
linear  case  [4],  the  fundamental  difference  between  the  long-  and  the  short-crested  sea  model  is 
that  a  wave  field  rather  than  a  wave  profile  must  be  dealt  with.  The  result  is  that  the  linear 
and  quadratic  response  functions  which  define  the  ship  dynamics  become  functions  of  vectors 
rather  than  scalars.  Given  a  directional  wave  spectrum  and  the  linear  and  quadratic  response 
functions  corresponding  to  some  response  of  interest,  the  model  provides  estimating  formulae 
for  a  variety  of  statistics.  When  the  nonlinearities  are  negligible,  the  estimates  reduce  to  those 
long  known  [4].  Relative  to  the  corresponding  formulae  for  the  long-crested  quadratic  model, 
there  are  no  real  surprises  in  the  short-crested  estimation  formulae,  only  the  greater 
complication  that  would  be  expected  because  of  the  increase  in  the  dimensions  of  the  problem. 
It  was  found  [3]  that  in  almost  all  cases  there  is  no  necessity  to  work  in  the  encounter  domain. 
This  result  implies  that  the  previous  restriction  of  the  long-crested  case  to  the  encounter 
frequency  domain  may  be  largely  removed.  The  formula  for  the  mean  value  of  the  response  in 
short-crested  seas  is  especially  simple,  in  that  it  was  found  to  be  hardly  more  computationally 
demanding  than  estimating  the  variance  of  a  purely  linear  system  in  short-crested  seas. 

Hasselmann  [1]  remarks  that,  since  the  basic  model  is  a  truncated  functional  series,  for 
mathematical  consistency  the  excitation  must  be  allowed  the  same  degree  of  nonlinearity  as  the 
system  being  excited,  and  thus  the  model  of  Dalzell  [3]  is  inconsistent  in  this  regard.  Thus  the 
overall  objective  of  the  present  investigation  was  to  remove  the  inconsistency  by  allowing  the 
wave  field  to  have  quadratic  nonlinearities.  The  detailed  present  objectives  included 
clarification  of  how  one  would  actually  carry  out  a  practical  computation,  what  further 
assumptions  are  being  made,  and  perhaps  most  importantly,  whether  the  complication  affects 
estimates  of  the  mean  nonlinear  response  which  is  the  most  important  aspect  of  quadratic 
response  in  practical  applications. 

The  inconsistency  at  issue  involves  a  priori  ideas  about  the  magnitude  of  quadratic 
nonlinearities.  In  most  if  not  all  hydromechanical  analyses  of  response  the  quadratic 
nonlinearities  are  “second  order”  effects.  Similarly,  quadratic  nonlinearities  in  the  wave  field 
have  long  been  known  to  be  “second  order”  from  the  Stokes  expansion.  The  question  is 
whether  “small”  nonlinearities  in  wave  field  have  a  disproportionate  influence  on  “small” 
nonlinearities  of  response.  For  engineering  purposes  the  answer  will  depend  on  the  details  of 
the  problem  to  a  greater  or  lesser  extent.  Accordingly,  a  mathematically  consistent  approach  is 
desirable,  even  if  it  merely  yields  indirect  estimates  of  the  magnitude  of  the  error  resulting 
from  the  inconsistency. 

Hasselmann  [1]  suggested  an  approach  to  the  consistency  problem  which  is  grounded  in 
the  Stokes  expansion,  dates  at  least  to  Tick  (5],  and  was  used  by  Hasselmann  et  al.  [6]  in  the 
interpretation  of  wave  bispectra.  More  recently,  the  approach  has  been  used  for  the  long- 
crested  case  by  Hineno  [7].  This  approach  centers  on  the  assumption  that  there  exists  a 
stationary  and  homogeneous  zero  mean  Gaussian  and  linear  “wave  field”  which  generates  the 
true  wave  field,  correct  to  second  order,  through  quadratic  wave-wave  interactions.  The  true 
wave  field  in  turn  then  excites  the  ship  system.  In  the  linear  context  the  approach  is 
conceptually  equivalent  to  the  often  used  model  in  which  the  response  is  visualized  as  the  result 
of  cascading  white  noise  through  a  wave  spectral  filter  and  exciting  the  ship  system  with  the 
result. 


QUADRATIC  WAVE-WAVE  INTERACTIONS 

The  solution  to  the  quadratic  wave-wave  interaction  problem  via  the  Stokes  expansion  has 
been  available  for  some  time,  having  been  used  by  Tick  [S],  Phillips  [8],  Longuet-Higgins  [9], 
and  Hasselmann  [6],  [10]  among  others.  Appendix  A  summarizes  the  results  in  notation  close 
to  that  used  by  Longuet-Higgins. 

For  present  purposes  the  important  aspect  of  the  hydrodynamic  analysis  is  that  the 
corrections  to  an  assumed  first  order  wave  field  due  to  second  order  wave-wave  interactions 
assume  the  form  of  the  quadratic  part  of  the  response  of  the  linear  plus  quadratic  functional 
model.  To  be  specific,  4L(X,t)  will  denote  a  deterministic  wave  field  composed  of  N 
superimposed  small  amplitude,  first  order  progressive  gravity  waves,  where  X  is  a  vector 
position  relative  to  an  X-Y  coordinate  system  fixed  in  the  mean  water  plane,  and  t  is  time. 

This  same  representation  of  the  wave  field  was  used  previously  [3]  and  may  be  written 

H 

fL(X,t)  =  ^tfeJbjExptiiCj.X-iWjt]}  (1) 

i 

where  bj  =  ajExp[Uj]. 

In  this  representation  aj  and  cj  are  the  amplitudes  and  phases  of  the  discrete  components  of  the 
small  amplitude  waves,  xj  is  the  corresponding  vector  wave  number,  and  coj  denotes  wave 
frequency.  The  wave  frequency,  coj,  is  assumed  to  be  connected  with  |ij|  via  a  dispersion 
relation  which  will  be  taken  to  be  the  deep  water  relationship  for  present  purposes. 

To  second  order,  wave-wave  interactions  will  distort  the  assumed  field  of  Eq.  1,  and  to 
this  order  the  corrected  or  “true”  wave  field  may  be  written: 

H 

?T(X,t)  -  ^RajbjExpCiZj^Z-iMjt]} 

J 

N  IT 

+ 2  b* 

J  k  _  , 

Exp[i(«j  +  Kk)  •  X  -  i(»j  +  o>k)t]f 

If  H 

+  —  Rf(/Cj,/ck) 

j  k  _  > 

Exp[i(/Cj  -  «k)  •  X  -  i(«j  -  (2) 

The  detail  of  this  transition  is  contained  in  Appendix  A,  and  explicit  formulae  for  the 
functions  and  R^*j,xk)  are  given  there  for  the  deep  water  case.  Except  that  they  are 

purely  real,  these  two  functions  have  the  same  symmetries  as  the  quadratic  frequency  response 
functions  previously  defined  [3],  which  relate  the  quadratic  sum  and  difference  frequency 
response  of  a  ship  system  to  the  amplitudes  of  component  small  amplitude  waves  in  the 
excitation.  In  this  case  the  functions  relate  the  components  of  the  true  field,  which  correspond 
to  sums  and  differences  of  vector  wave  numbers,  to  the  interactions  of  the  constituent  small 
amplitude  waves. 

The  analogy  to  the  ship  system  response  case  would  be  much  closer  if  the  first  term  of 
Eq.  2  were  different  from  Eq.  1  by  a  factor  of  some  function  of  wave  number  inside  the 
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summation.  Such  a  factor  would  correspond  to  a  linear  frequency  response  function.  In  this 
case  the  hydrodynamic  theory  shows  that  there  is  no  linear  transformation  of  the  small 
amplitude  wave  field,  and  thus  the  analogy  must  be  completed  by  taking  the  linear  frequency 
response  function  to  be  unity. 

Although  it  has  not  been  done  explicitly  in  Appendix  A,  it  appears  that  the  potential  of 
the  true  wave  Held  may  also  be  expressed  in  a  form  similar  to  that  of  Eq.  ?.  This  result  tends 
to  justify  the  idea  that  the  properties  of  the  true  wave  field  to  second  order  may  be  represented 
as  the  response  of  a  physical  linear  plus  quadratic  system  to  a  Active  “linear”  small  amplitude 
field.  It  also  strongly  suggests  that  the  theoretical  response  of  the  ship  system  to  the  true  field 
may  be  approached  as  the  result  of  one  linear  plus  quadratic  system  (the  ship)  acting  upon 
another  (the  wave-wave  interactions). 

REPRESENTATION  OF  THE  WAVE  FIELD 

The  hydrodynamic  arguments  suggest  that  the  true  wave  field  may  be  taken  to  be  a 
functional  of  a  “linear”,  small  amplitude  field,  at  least  to  second  order.  In  accordance  with 
the  previous  work  involving  system  response  [3],  if  the  linear  field  is  continuous  and 
differentiable  and  vanishes  outside  some  finite  region  in  the  space  of  all  admissible  linear 
fields,  or  if  the  field  is  square  integrable,  the  true  wave  field  may  be  represented  as  a 
functional  Taylor  expansion  in  the  neighborhood  of  zero  wave  field.  With  further  general 
assumptions  noted  in  Appendix  B,  the  Taylor  expansion  through  the  second  order  term  may  be 
written  as  follows: 

CT(t,X,Y)  =  f  L(t,X,Y) 

+/-/« * 

fL(t-ti  .X-Xi  ,Y— y j) f  L(t— ta  ,X— xa  ,Y— ya)  x 

dtj  dxj  dyt  dt8  dxg  dyg  (3) 

where  the  integrals  are  infinite.  The  kernel  within  the  integral  in  Eq.  3  is  a  function  of  time 
and  space  differences  and  contains  the  physical  relationship  between  the  linear  field  and  the 
second  order  components  of  the  true  wave  field.  The  variables  t,  X,  and  Y  denote  absolute 
time  and  position  in  a  fixed  frame  of  reference.  As  in  the  previous  work  it  has  been  assumed 
that  the  integrands  differ  from  zero  only  in  a  finite  range  of  the  dummy  arguments. 

A  difference  between  Eq.  3  and  the  corresponding  representation  [3]  of  the  response 
should  be  noted.  The  response  of  the  ship  system  is  supposed  invariant  with  absolute  position, 
and  in  the  representation  of  the  response  this  restriction  was  taken  into  account  in  the 
formulation  for  the  encounter  domain.  In  effect,  had  the  response  been  presumed  a  function 
of  absolute  position,  the  form  of  the  final  functional  representation  would  have  been  similar  to 
Eq.  3  except  that  a  linear  convolution  would  also  have  appeared. 

The  hydrodynamics  of  the  problem  noted  in  the  last  section  shows  that  there  is  no  linear 
transformation  of  the  linear  field.  This  fact  accounts  for  the  absence  of  a  linear  convolution  in 
Eq.  3. 
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An  immediate  problem  which  arises  with  this  representation  is  both  theoretical  and 
practical  (from  a  computational  point  of  view).  The  quadratic  wave-wave  interaction  response 
functions  are  known  in  the  wave  number  domain.  The  consequence  of  the  representation  of 
Eq.  3  is  that  the  second  degree  kernel  and  the  quadratic  response  functions  must  form  special 
multi-dimensional  Fourier  transform  pairs.  Unfortunately  the  analysis  of  Appendix  A  shows 
that  the  quadratic  response  functions  do  not  possess  all  the  required  properties  [11]  to  make 
the  required  transforms  valid.  The  problem  with  these  functions  is  that  they  are  unlimited  as 
wave  number  tends  to  infinity. 

A  way  out  of  the  problem  is  detailed  in  Appendix  B.  Essentially,  a  wave  number  window 
is  hypothesized  such  that  the  final  expansion  represents  the  interactions  only  between  wave 
components  having  wave  numbers  less  than  some  large  but  finite  value.  The  window  is 
represented  by  a  linear  transformation,  and  the  revised  formulation  may  be  summarized  in  a 
cascade  diagram: 


fL(t.x) 

gt(.) 

CT(t,X) 

fl(.) 

g8() 

Here  the  notation  £N(t,X)  represents  the  nominal  linear  wave  field.  This  field  is  presumed  to  be 
operated  upon  by  a  linear  system  fj(.)  in  the  time/space  domain  to  produce  a  linear  small 
amplitude  wave  field  £L(t,X)  which  is  wave  number  limited  but  defined  as  before.  The  linear 
small  amplitude  field  is  acted  upon  by  the  quadratic  wave  interaction  system  with  kernels  g,(.) 
and  g2(.)  to  result  in  the  true  wave  field  £j<t,X).  The  fast  functional  relation  is  the  same  as 
Eq.  3,  with  the  missing  linear  kernel  presumed  to  be  a  delta  function. 

The  Taylor  expansion  of  the  true  field  in  terms  of  the  nominal  field  becomes 

fT(t,X,Y)  =  fff  h1(t1.x1,y1)fs(t-t1,X-x1,Y-y1)dt1dx1dy1 

+/.../h8(t1,x1,y1,t8,x8.y8)  x 

CN(t-t1  ,X-Xl , Y-y  f  M(t— 18  ,X-x8  ,Y-y 2)  x 

dtj  dX|  dyx  dt8  dx8  dyg  (4) 

and  the  response  functions  relating  nominal  and  true  wave  fields  may  be  written  as  follows: 


Tt0c)  =  /../hjCtpXj.yj)  Expf-KwxXi+KTy^oti)]  dttdXldy 

=  /../hs(t1,x1,y1,tB,x8,y8)  x 

ExpC-U/Cxx^+WTiy^o^tJ]  x 

Exp[-i(A:XaX3+A:T8y8-a>3t8)]  dt^.-dyg 


l 


(5) 


(6) 
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rt(*i ,*b)  =  /••/hB(ti.x1,y1,te,x8,y8)  x 

ExpE-^/Cx^i+CTiy^W^j)]  x 


Exp[i(/cXBx8+/cT8y8-w8ta)]  dt1...dya  (7) 

where  frequency  has  been  omitted  from  the  arguments  of  the  response  functions  because  wave 
frequency  is  a  function  of  wave  number  via  the  dispersion  relation. 

From  the  cascade  analysis  of  Appendix  B  the  relationships  between  the  functions  defined 
by  Eqs.  5  through  7,  the  quadratic  wave-wave  interaction  functions  defined  previously,  and  the 
linear  window  response  function  TN(x)  are  defined  as  follows: 

Tt(«)  =  T„(«)  (8) 

=  7n(*i)  (9) 

®T(|Ci>*a)  =  Tif(/Ci)  Tk(k8)  (10) 

The  relation  between  the  window  response  function  and  the  space/time  domain  window 
fj(.)  is  the  same  as  Eq.  5  with  fj(.)  substituted  for 

Except  that  it  must  approach  zero  at  a  sufficient  rate  to  make  the  left  hand  sides  of  Eqs. 

9  and  10  approach  zero  for  large  wave  number,  the  window  response  function  is  quite 
arbitrary.  In  what  follows  it  will  be  assumed  to  be  analogous  to  a  non-realizable  low  pass  filter 
with  unity  gain  up  to  some  cutoff  wave  number  and  a  rapid  transition  to  zero  above.  Since  the 
wave-wave  interaction  functions  are  real,  this  assumption  makes  Eqs.  8  through  10  purely  real. 

The  interpretation  of  the  windowed  response  functions  just  defined  is  identical  to  that  of 
the  last  section.  If  the  nominal  wave  field  is  represented  by  a  sum  of  cosinusoids  exactly  as  in 
Eq.  1  and  this  representation  is  substituted  into  Eq.  4,  the  resulting  expression  for  the  steady 
state  true  wave  field  becomes 

N 

?T(t,X)  =  ^  Tt(7cj)  Re\ b  j  Exp[isj  ♦  X  -  icujt  ]| 
i 

N  » 

+  ^  QT(/Cj  »**)  * 

J  k  <  -  , 
i?e|bj  bk  ExptU/ij  +  «k)  ♦  X  -  i(co j  +  a>k)t]j 

N  N 

+  x 
t  X  c  * 

fle|bjbk  Exp[i(«j-/ck)»X-i(uJ-«k)t]j  (11) 

As  before,  the  windowed  quadratic  response  functions  Qj(icj,jek)  and  RT(5cj,jck)  relate  those 
components  of  the  nominal  field  which  correspond  to  sums  and  differences  of  vector  wave 
numbers,  to  the  interactions  of  small  amplitude  waves  which  reside  within  the  window.  The 
range  of  wave  numbers  present  in  the  linear  part  of  the  true  field  is  determined  by  Tx(xj)  in 
the  first  term  of  the  expression. 
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SPECTRAL  REPRESENTATION  OF  THE  TRUE  WAVE  FIELD 

The  next  step  in  the  development  involves  the  introduction  of  randomness.  This  step  may 
be  approached  in  the  same  manner  as  in  the  case  of  response  to  a  linear  wave  field  [3].  For 
this  purpose,  it  is  assumed  that  the  nominal  linear  wave  field  £N(X,Y,t)  is  dispersive, 
homogeneous  in  space,  stationary  in  time,  and  Gaussian.  It  is  further  assumed  that  the 
elevation  variance  of  this  nominal  wave  field  may  be  represented  by 


Var(f„)  =  ff  FKUX  ,*T)d*x  d/cT  (12) 

where  Fn<xx»xy)  is  a  vector  wave  number  spectrum. 

These  assumptions  together  imply  that  the  field  has  a  “spectral  representation”  [12],  [13], 
[1].  Although  this  representation  is  often  written  in  the  form  of  a  Fourier-Stieltges  integral,  it 
is  convenient  to  work  with  the  corresponding  approximating. Fourier  sum.  Accordingly,  the 
approximating  Fourier  sum  of  the  spectral  representation  of  the  nominal  field  is  taken  to  be 

f„(t,X)  =  2j?e{w(j)Exp[i*XJX  +  i/cTjY-i«jt]j  (13) 

J 

where  W(j)  is  a  complex  random  function  of  the  interval  Axxj,AxYj.  The  initial  assumptions 
that  the  nominal  field  is  stationary  and  homogeneous  and  that  the  process  is  Gaussian  imply 
that  the  random  functions  are  Gaussian,  zero  mean,  and  statistically  independent  for  different 
wave  numbers  and  directions  of  propagation.  Under  these  assumptions  the  ensemble  expected 
values  of  the  function  and  its  various  order  moments  may  be  defined  and  become  either  zero 
or  relatively  simple  functions  of  the  variance  spectrum.  Appendix  C  summarizes  the  results 
through  sixth  order  moments. 

As  in  the  deterministic  case  of  the  previous  sections,  a  substitution  of  Eq.  13  into  Eq.  4 
yields  the  statistically  steady  state  spectral  representation  of  the  true  wave  field  in  the  form  of 
an  approximating  Fourier  sum.  Because  the  fandom  functions  are  analogous  to  the  complex 
amplitudes  of  Eq.  11,  the  algebra  of  the  substitution  and  subsequent  reduction  is  identical  to 
that  of  the  deterministic  case.  The  resulting  spectral  representation  of  the  true  wave  field  is 
identical  in  form  to  Eq.  1 1  and  may  be  written 

fT(t,X)  =  ^  Tx(7Cj) /?ejw(j)  ExptiVX-ic^Oj 
J 

+  i'2  2  x 

*e|w(j)W(k)  Exp[i(«j  +  •  X  -  Kwj  +  &>k)t]} 

+  iS  x 

*ejw(j)W*(k)  Exp[iOcj-/ck)*  X  -  -  &>k)t]j  (14) 
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where  it  should  be  noted  that,  as  in  all  the  development  thus  far,  the  spatial  coordinate  system 
is  fixed,  and  the  frequencies  are  wave  frequencies  related  to  scalar  wave  number  through  a 
dispersion  relation. 

STATISTICAL  PROPERTIES  OF  THE  TRUE  WAVE  FIELD 

The  true  and  nominal  wave  fields  differ  by  the  windowing  function  and  the  effects  of 
wave-wave  interaction.  The  spectral  representation,  Eq.  14,  and  the  properties  of  the  complex 
random  functions  noted  in  Appendix  C  allow  some  of  the  important  statistical  properties  of 
the  true  field  to  be  inferred. 

First  among  the  important  statistics  is  the  statistical  mean  value.  When  the  ensemble  mean 
value  of  Eq.  14  is  taken,  almost  all  the  contributing  terms  drop  out  by  virtue  of  Eqs.  C.3,  C.4, 
and  C.6.  The  value  of  the  single  term  remaining  is  determined  by  Eq.  C.S,  and  the  resulting 
mean  value  of  the  true  wave  field  may  be  written 

(Ct(I.X))  —  ]  Rt(*j  »*])  Fw(<Cjrj  ,*Tj)  AcXJ  AacTj  =  0  (15) 

i 

where  the  notation  <..>  denotes  the  ensemble  expected  value,  and  the  expression  is  zero  Finally 
because  the  wave-wave  interaction  component  of  RxOcj.xj)  is  zero  (Appendix  A).  Thus  Eq.  15 
expresses  the  physically  expected  result  that  the  true  wave  field  is  a  zero-mean  process. 

The  other  statistics  of  primary  interest  are  the  true  wave  variance  and  the  true  wave 
number  spectrum.  Roth  may  be  approached  by  first  developing  the  space-time  auto  covariance 
function  <£x(t,X)  {jft  +  t,X + p)>  in  which  r  denotes  a  time  difference  and  q  a  vector  position 
difference.  Using  the  spectral  representation,  Eq.  14,  for  each  of  the  factors  inside  the 
expectation  operator,  expanding  the  real  part  operation  into  a  sum  of  the  function  and  its 
complex  conjugate,  and  forming  the  required  product  results  in  thirty-six  terms  involving 
multiple  summations  of  complex  exponentials,  the  response  functions,  and  the  complex  random 
amplitudes.  Once  the  expectation  is  taken,  all  but  twelve  of  these  products  drop  out  by  virtue 
of  Eqs.  C.7  and  C.8,  and  the  remainder  are  evaluated  and  simplified  by  Eqs.  C.4,  C.S,  C.9, 
and  C.10.  After  a  further  algebraic  reduction,  which  parallels  that  described  for  the  response 
[3],  the  approximating  Fourier  sum  representation  of  the  space-time  auto  covariance  of  the  true 
wave  field  becomes 

(CT(t,X)  fx(t+T,X+p))  = 

Z  Tt(*j)fk(*XJ  »*Tj)  Exp[iKj  •  p-iUjT  ]}  A*XJ  A/Cyj 
J 

+Zj  Xy  **k)  fn(*xj  •*  Tj)Flt(*Xk  »*Tk)  x 

/?e|Exp[i(c j  +  *k)  •  p  -  l(ci) j  +  &>k)r]|  A/cXJ  A«T j  A#cXk  A/cTk 


J  k 


Rt(*j  •*!*)  Fk(*XJ  t*Tj)F lt(*Xk  »*Tk)  x 
RejExp[i(JCj  -  #ck)  •  p  -  i(a»j  -  «k>T]j  A/cXJ  A*TJ 


A*Xk 


A/cTk 
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At  this  stage  all  the  complex  random  amplitudes  have  been  eliminated  and  with  them  the 
necessity  for  the  Fourier  sum  approximations  to  Fourier-Stieltges  integrals.  Thus  when  the 
range  of  the  summation  is  allowed  to  increase  and  the  wave  number  intervals  to  decrease,  the 
limiting  form  involves  ordinary  integrals 

(CT(t,X)fT(t+r,X+p))  =  /../t?0c)  Fm0c)  Re\  Exp[i*  •  p  -  i«r  ]}  die 

fle jExpfU*!  +  *B)  •  p  -  i(«l  +  &>8)t]|  d«i  dicg 

+/../R!(<Ci,«a)PM(«1)P1(0cg)  x 

Re jExptK^i  -  icB)  •  p  -  !(&>!  -  a>B)r]|  d*t  dicB  (17) 

Now  if  the  time  and  vector  space  differences  are  set  to  zero,  Eq.  17  represents  the  true 
wave  variance 

Var(fT)  =  /../T|0c)PM(3c)dic 

+/"/iQ t<*i  •*«)  +  rt (*i  .*e){  (18) 

Since  the  response  functions  defined  by  Eqs.  8  through  10  are  here  assumed  to  be  real  valued, 
this  result  is  exactly  that  which  would  be  expected  from  the  previous  development  [3]. 

In  accordance  with  convention  (8]  the  vector  wave  number  spectrum  of  the  true  wave  field 
FT(ic)  may  be  defined  as  follows: 

(?T(X)  fT(X+p))  =  y*..y*FT(7c)Exp[i**  p]d*  (19) 

where  the  wave  field  is  assumed  to  be  homogeneous.  An  expression  for  the  wave  number 
spectrum  of  the  true  wave  field  in  terms  of  that  of  the  nominal  wave  field  is  obtained  by 
setting  the  time  difference  in  Eq.  17  to  zero  and  making  changes  of  variable  so  as  to  put  the 
result  in  the  form  of  Eq.  19.  By  this  process  the  expression  for  the  true  wave  number  spectrum 
becomes 

Ft(*)  =  T|(«)  F„(*) 

+  f  Q?[(*+v)/2,(«-v)/S]  F„[(*  +  v)/2]  F„[(k-v)/3]  dv 

+-^-/..y*R|[(«+v)/8  ,(v-«)/3]  Fk[(k  +  v)/2]  F„[(v-«)/8]  dv  (20) 

The  convolutions  in  this  expression  are  of  the  general  form  which  would  be  expected  from 
the  time-only  theory  for  linear  plus  quadratic  systems.  In  the  practical  case  one  might  want  to 
analyze  an  observed  wave  number  spectrum  by  deriving  the  nominal  spectrum.  The  results  just 
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shown  illustrate  the  difficulties  in  such  an  endeavor  —  the  nominal  vector  wave  number 
spectrum  is  defined  by  integral  equation  of  greater  complexity  than  that  which  obtains  in  the 
time-only  case.  The  present  model  is  oriented  more  toward  synthesis  than  analysis,  although 
some  analysis  steps  have  been  made  [<$]. 

Because  in  marine  practice  the  directional  variance  spectra  of  waves  are  usually  considered 
in  the  wave  frequency/direction  form,  it  was  of  interest  to  see  whether  an  equivalent  or  simpler 
form  for  the  true  wave  spectrum  would  be  obtained  in  the  wave  frequency/direction  domain. 
The  required  transformations  were  carried  through  far  enough  to  see  that  the  equivalent 
expression  would  be  far  more  complicated.  Accordingly,  if  the  computations  equivalent  to  Eqs. 
18  or  20  need  to  be  done  in  this  domain,  it  appears  less  troublesome  to  transform  the  nominal 
frequency/direction  spectrum  to  the  vector  wave  number  domain,  carry  out  the  convolutions, 
and  then  transform  the  true  vector  wave  number  spectrum  back  to  the  wave 
frequency/direction  domain. 

REPRESENTATION  OF  THE  RESPONSE 

The  representation  of  the  ship  response  to  the  true  wave  field  is  essentially  the  same  as  in 
the  previous  work  (3).  It  is  assumed  that  the  ship  is  proceeding  at  mean  speed  U  along  the  X 
axis  in  a  fixed  X-Y  frame  of  reference.  Component  wave  directions  will  be  denoted  by  ji.  and 
the  resulting  conventions  are  indicated  in  the  figure. 


Ail  (x-y)  coordinate  system  moving  with  the  ship  is  conventionally  defined  as  follows: 


x  -  X-Ut 

y  -  Y 

and  accordingly  the  encounter  frequency  corresponding  to  a  given  wave  component  becomes 
o  *■  o>  -  *XU 

—  encounter  frequency 

As  in  the  previous  work,  the  ship  response  |(t)  may  be  represented  by  a  functional  Taylor 
expansion  in  the  neighborhood  of  zero  true  wave  field  as  encountered  by  the  ship.  The 
resulting  model  is  written  as  follows: 

t(0  *  ///k1(ti.*i»y1)CT(t-ti,x1.y1)dt1dx1dy| 

+/••  /*,( t1  ,xx  ,y  t  ,t8  ,x8  ,yg)  x 

f T(t-t4  ,xx  ,y j)  f T(t-t,  ,x8  ,y8)  dt4  dXj  dyt  dt8  dx,  dyg  (21) 

The  first  (constant)  term  in  the  expansion  has  been  dropped  on  the  presumption  that  no 
response  of  current  interest  occurs  in  the  absence  of  waves,  and  the  series  has  been  truncated 
at  the  quadratic  term.  The  convolutions  are  with  respect  to  time  and  space  differences.  Since 
this  is  an  encounter  domain  model,  the  space  differences  are  with  respect  to  a  reference  point 
on  the  moving  ship.  The  linear  ship  dynamics  are  assumed  to  be  represented  by  the  linear 
kernel  k((t|,x1,yi)  and  the  quadratic  nonlinearities  by  the  kernel  k2(t],Xj,y|,t2,X2,y2)-  Apart 
from  the  explicit  note  that  the  local  wave  field  is  the  true  Held,  Eq.  21  is  identical  to  Eq.  3  of 
the  previous  work  [3]. 

The  identification  of  linear  and  quadratic  response  functions  corresponding  to  the  kernels 
of  Eq.  21  is  identical  to  that  of  the  previous  work.  The  steady  state  response  to  assumed 
“true”  wave  fields  composed  of  one  and  two  superimposed  linear  regular  waves  is  obtained, 
and  in  this  process  the  relationships  between  the  kernels  and  response  functions  may  be  defined 
in  the  following  form: 

Tt(*,a)  =  f..f  k^tj.xjjj)  ExpOKKxXi+w^  +  ot!)]  dt1dx1dyl  (22) 

Q^j.ffi.ca.o,)  =  y*”#/*kB(t1,x1,y|,t#,x8,yg)  x 

Exp[-i(*xlx,+*T1yl  +  o1tl)]  x 

Exp[-i(*x8x8  +  ieT8yg+o8t8)]  dtt...dys  (23) 
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*  /../k.(t1.xl,y1.t„x..y.)  x 

ExpC-KiCxiXt+cr^^ait^]  x 

Exp[i(xXIx*+*Ttyg+a1tt)]  dt^-.dy,  (24) 

As  before  [3],  Eq.  22  defines  the  linear  response  function,  Eq.  23  the  quadratic  sum 
frequency  response  function,  and  Eq.  24  the  quadratic  difference  frequency  response  function. 
The  existence  of  these  transforms  requires  that  the  kernels  be  square  integrable  —  in  effect  that 
the  time/ space  “memory”  of  the  physical  ship  system  be  finite  (or  that  some  artifice  to  make 
this  so  can  be  introduced  if  necessary).  The  notation  departs  from  that  used  previously  (3)  in 
that  the  encounter  frequency  corresponding  to  each  vector  wave  number  argument  is  noted 
explicitly  in  the  argument  of  each  function.  In  the  previous  work  if  a  dispersion  relation  and 
the  ship  speed  are  assumed,  the  encounter  frequency  is  determined,  and  accordingly  the 
encounter  frequency  argument  was  dropped  on  this  basis.  In  portions  of  the  present  work  this 
understanding  is  not  necessarily  valid,  and  the  present  more  cumbersome  notation  is  of 
advantage.  However,  the  meaning  of  the  response  functions  is  identical  to  that  previously 
given.  The  linear  response  function  is  the  normalized  amplitude  of  response  in  the  encounter 
frequency  to  a  linear  regular  wave  component.  The  sum  (difference)  frequency  response 
function  is  the  normalized  amplitude  of  response  at  an  encounter  frequency  equal  to  the  sum 
(difference)  of  the  encounter  frequencies  corresponding  to  each  of  two  linear  regular  wave 
components. 

From  the  definitions  two  symmetry  relations  arise: 

Q|<*i."i.*a.®a)  =  Q^**.**.*!.*!) 

Ht(«c1,cr1,#cs.ara)  =  Rt  v)  (25) 

as  well  as  the  formal  relation 

(26) 

SPECTRAL  REPRESENTATION  OF  THE  RESPONSE 

Because  the  response  model,  Eq.  21,  is  for  the  encounter  domain,  and  the  spectral 
representation  of  the  true  field,  Eq.  14,  involves  the  fixed  spatial  coordinates  of  Figure  1,  a 
conversion  is  necessary.  Transformation  of  Eq.  14  from  the  fixed  (X-Y)  coordinates  to  the 
moving  (x-y)  system  results  in  the  following  spectral  representation  of  the  true  wave  field: 

M4'*)  "  ^  Tt02,)  *•}*())  Exp[i«rX-iojt]j 
J 

+  x 

ff*jw(J)W(k)Exp[i(3cJ  +  Ek).X-i(oJ  +  ok)t]j 
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(27) 


J  *#|w(j)W*(k)  Bzp[i(Xj  —  Sk)  •  X  —  i(<Tj  —  ok)* 

The  form  of  Eqs.  14  and  27  is  the  same.  Wave  frequencies  in  Eq.  14  are  replaced  by  encounter 
frequencies,  and  the  fixed  position  vectors  become  relative  to  the  mean  ship  position.  The  only 
difference  between  this  stage  and  the  corresponding  stage  in  the  previous  work  [3)  is  the 
presence  of  wave-wave  interaction  terms.  Equation  27  is  the  Fourier  sum  approximating  the 
local  random  wave  field  as  it  would  be  sensed  aboard  the  ship. 

The  interest  here  is  in  the  statistically  steady  state  ship  response  to  the  true  wave  field.  As 
in  the  previous  work,  the  substitution  of  the  Fourier  sum  spectral  representation  of  the  true 
wave  field,  Eq.  27,  into  the  time  domain  model,  Eq.  21,  yields  a  spectral  representation  of  the 
steady  state  response  in  the  form  of  an  approximating  Fourier  sum.  It  was  anticipated  from 
the  form  of  the  two  equations  that  this  direct  substitution  would  yield  a  result  which  in  part 
would  involve  contributions  corresponding  to  nonknearities  of  higher  degree  than  those  which 
are  of  interest  here.  The  approach  taken  was  to  ignore  the  potential  inconsistency  pending  the 
complete  result. 

The  mechanics  of  this  substitution  involves  replacing  the  real  pans  of  the  various 
expressions  in  Eq.  27  by  half  the  sum  of  the  respective  expression  and  its  complex  conjugate, 
substituting  the  result  in  Eq.  21,  and  carrying  out  the  integrations  with  respect  to  the  time  and 
space  difference  variables  by  applying  the  response  function  definitions,  Eqs.  22  through  24. 
The  raw  result  involves  6  summations  from  the  first  integral  of  Eq.  21  and  36  summations 
from  the  second.  These  summations  occur  in  conjugate  pairs  so  that  the  final  results  may  be 
expressed  in  real  part  notation.  The  form  of  the  summations  allows  some  terms  to  be 
combined,  and  the  symmetry  relations,  Eqs.  25  and  26  serve  to  combine  others.  The  net  result 
is  12  terms  of  different  form  which  may  be  expressed  as  8  summations,  and  the  final  result  is 
given  as  Eq.  28  on  the  next  page. 

It  will  be  remembered  that  the  function  T^j)  is  a  wave  number  window  which  has  been 
assumed  to  be  purely  real  and  whose  purpose  is  to  make  the  wave-wave  interaction  model 
valid.  This  window  function  has  been  combined  with  the  wave-wave  interaction  response 
functions  and  R^Xj.Xn,)  in  Eqs.  9  and  10.  The  cascaded  response  functions  defined 

by  Eqs.  9  and  10  have  been  used  in  the  immediately  previous  development  to  compress 
notation.  In  the  case  of  Eq.  28  the  window  functions  have  been  factored  outside  of  the  real 
part  operator.  All  the  summations  of  Eq.  28  are  over  wave  number,  and  the  net  result  is  as 
would  be  expected;  that  is,  in  any  of  the  summations,  any  term  which  involves  a  wave  number 
above  the  cutoff  determined  by  the.  window  function  will  be  attenuated  or  zeroed  entirely. 

Each  of  the  summations  of  Eq.  28  has  a  similar  form:  first,  the  product  of  as  many 
window  functions  as  there  are  summations,  and  then  the  real  part  of  a  complex  expression 
made  up  of  factors  of  three  kinds.  Within  the  complex  expression  in  each  case  is  the  product 
of  as  many  complex  random  amplitudes  W(j)  as  there  are  summations.  These  factors  are 
followed  by  an  exponential  factor  containing  time  and  sums  and  differences  of  encounter 
frequencies.  Finally,  the  last  factor  involves  products  and  sums  of  wave-wave  interaction  and 
system  response  functions. 
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t(t)  -  ^T1(*J)*.jW<J)Ixp[iajt]Tt(*J,«rJ)j 

♦tT!  S  T.(*j)  T*(Em)  W*(m)  Bxp[i(«7J+<r-)t]  x 

[Qf  T^Xj+TC.j.aj+a.)  +Qe(*j,<r  ]j 

H££T"(*')T*(*“)**^J)W(m)Exp[i(<FJ‘<r«)t)  X 

Tt(*j-i5m,aj-am)  +Rt(*J,aJ,*m,<rm)  ]\ 


+£]£  £  £ T*(ZJ) T«<*.>*#|wKj) "*(**) **U)  * 

Exp[i(ffJ+am+an)t]Qf(*J.*m)Qt(*J+*BI,<r,+am,*m,«rn)]j 

+£EEEt  «<*,)  TM(7Em)  !•,(**)  J?#{w*(J)  W*(m)  W(n)  x 

Ixp[i(<rJ+am-aB)t][Qc(*J.*m)R|(XJ+*m,aJ+am,7EB.alk) 

+Rt(Ej,E0)  Qt(«J-**,aJ-«rB,Ein.am) 


+i£  E  E  E  T«<*»>  *><*->  *.<*•>  T-(*»> x 

jw*(J)  W*i(m)  W*l(n)  W*(p)  Exp[i(<Tj+am+aB+a p)t]  x 

+i£  £  £  £  t»(2j)  T»<  v * 

*®|**(J)  W*(m)  W*<n)  W(p)  Exp[l(a,+ffm+aB-ff,)t]  x 
Qf  (2j.2m)  Rf  (*».*,)  Qt(*j+*m.*  j+»m>*n-*j,<yn-ap)  | 

+?EEEEt  »<*,) Tm^*)  Tp(*.) Th(Z,)  x 

W*(m)  W(n)  W(p)  Exp[i(aj+a1B-or0-ap)t]  x 
+  R((Kj,xB)R((Km.*p)Qt(ici-*B,aJ-*n,Jcm-Kp,<rm-(rp) 

+  Rf(*j.*p)  Rf(*n,2i»)  ]| 

(28) 
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The  degree  of  nonlinearity  of  each  term  of  the  expression  is  easily  determined  by  a  count 
of  the  number  of  summations.  In  each  term  the  number  of  wave  number  window  functions  in 
the  product  equals  the  number  of  summations.  The  number  of  complex  amplitude  functions  in 
the  product  is  the  same,  as  is  the  number  of  encounter  frequencies  summed  in  the  exponential. 
Had  the  Taylor  expansion  of  the  response,  Eq.  21,  been  truncated  after  the  fourth  instead  of 
after  the  second  functional,  the  corresponding  approximating  Fourier  sum  would  have  been  of 
very  similar  form  to  Eq.  28,  with  the  response  function  corresponding  to  the  N’th  degree 
kernel  appearing  in  an  N-fold  summation  which  also  involves  an  N-fold  summation  of 
encounter  frequencies. 

The  first  term  of  Eq.  28  is  the  only  one  of  first  degree;  it  represents  the  linear  ship  system 
response  to  the  nominal  wave  field  and  is  the  same  as  the  corresponding  term  of  the  response 
representation  of  [3]  except  for  the  window  function.  The  next  two  terms  are  second  degree 
and  represent  the  quadratic  response  in  sum  and  difference  encounter  frequencies.  Since  the 
final  factor  in  these  complex  expressions  is  a  sum  of  two  terms  of  the  form  Q{ficj,xm)» 

T$(*j  +  *m,°j  +  °m)  and  the  response  of  second  degree  has  two  qualitative 

origins.  One  is  the  quadratic  ship  system  response  to  pairs  of  linear  wave  field  components 
(identical  to  the  results  in  the  previous  work  [3]).  The  other  is  the  linear  ship  system  response 
to  all  the  wave-wave  interaction  sum  or  difference  frequency  components  present  in  the  true 
wave  field.  This  last  contribution  is  the  first  of  the  modifications  due  to  the  difference  between 
true  and  nominal  wave  fields. 

None  of  the  remaining  terms  in  Eq.  28  have  counterparts  in  the  earlier  work  [3].  The 
fourth  and  fifth  terms  in  Eq.  28  are  of  third  degree  and  involve  response  in  an  encounter 
frequency  equal  to  the  sum  as  well  as  a  partial  difference  of  three  linear  encounter  frequencies. 
These  terms  represent  all  the  possible  ship  system  quadratic  responses  to  all  the  pairs  of  true 
wave  field  components  which  may  be  formed  with  one  linear  wave  field  component  and  one  of 
the  components,  of  the  wave  field  produced  by  quadratic  wave-wave  interactions. 

The  sixth  through  eighth  summations  in  Eq.  28  correspond  to  fourth  degree  nonlinearities. 
This  part  of  the  response  occurs  at  sums  and  differences  of  four  linear  encounter  frequencies. 
These  terms  represent  all  the  possible  ship  system  quadratic  responses  to  pairs  of  the 
components  in  the  true  wave  field  which  are  produced  by  wave-wave  interactions. 

The  form  of  the  arguments  of  some  of  the  response  functions  in  Eq.  28  has  come  about 
purely  as  a  result  of  the  definitions,  Eqs.  22  through  24.  The  reason  for  the  addition  of  the 
encounter  frequency  argument  in  the  definition  of  these  functions  may  now  be  apparent.  For 
example,  in  the  linear  function  Tj(*j  +  xm, Oj  +  om)  t*ie  argument  *j  +  *in  defines  *  wave-wave 
interaction  component  which  may  have  a  propagation  direction  different  from  those  of  the  two 
linear  wave  field  components  involved  and  almost  certainly  a  different  wave  length.  The  wave 
frequency  which  would  be  inferred  by  applying  the  deep  water  dispersion  relationship  to  the 
absolute  value  of  the  vector  sum  is  not  necessarily  the  same  as  the  sum  of  the  wave  frequencies 
of  the  components,  nor  is  the  corresponding  encounter  frequency.  Since  many  ship  motion 
responses  involve  a  mixture  of  dynamic  resonances  with  a  geometric  matching  of  wave  shape 
and  vessel  shape  (wave  matching),  it  is  important  to  keep  track  of  such  differences. 

As  has  been  noted,  the  result,  Eq.  28,  has  been  obtained  under  assumptions  which  are 
mathematically  inconsistent  in  the  context  of  series  expansions.  In  particular,  it  has  been 
assumed  that  the  wave-wave  interactions  are  strictly  limited  to  the  quadratic  (“second  order”) 
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type;  that  is,  cubic  and  nonlinear  interactions  of  higher  degree  are  assumed  to  be  identically 
zero.  Similarly,  nonlinear  ship  system  interactions  of  degree  higher  than  two  are  assumed  not 
to  exist.  Equation  28  is  a  Fourier  sum  approximating  the  spectral  representation  of  the 
response  under  these  assumptions. 

It  is  helpful  to  define  an  “order” -(in  the  magnitude  sense)  of  the  various  terms  of  Eq.  28. 
The  most  convenient  way  is  to  assume  that  the  magnitude  of  response  is  most  closely 
associated  with  the  magnitude  of  the  wave  components.  If,  as  is  necessary  when  dealing  with 
weakly  nonlinear  systems,  it  is  assumed  that  the  components  of  the  nominal  (linear)  wave  field 
are  small,  then  products  of  components  are  smaller  yet.  In  the  present  context  the  complex 
amplitudes  W(j)  are  a  measure  of  the  magnitude  of  each  wave  component.  If  these  amplitudes 
are  small  and  assumed  to  be  of  “first  order”,  then  N-fold  products  of  the  complex  amplitudes 
are  of  N’th  order.  On  this  ordering  assumption  it  may  be  seen  that  the  order  of  the  various 
terms  in  Eq.  28  is  the  same  as  the  degree  of  nonlinearity.  Thus,  if  it  seems  permissible  to 
disregard  response  of  third  and  higher  order,  only  the  first  three  summations  of  Eq.  28  will 
remain. 


EXPECTED  VALUE  OF  EQUATION  28 

The  form  of  the  statistical  mean  value  of  Eq.  28  as  it  stands  is  instructive.  The  required 
algebra  is  not  overpowering,  and  thus  the  order  of  the  various  terms  will  be  temporarily 
ignored. 

The  ensemble  average  of  six  of  the  eight  summations  of  Eq.  28  may  be  immediately  seen 
to  be  zero  from  the  relations  of  Appendix  C.  In  particular,  the  mean  of  the  first  summation  is 
zero  by  Eq.  C.3,  and  the  mean  of  the  second  by  Eq.  C.6.  The  mean  of  the  fourth  and  fifth 
summations  is  zero  by  Eq.  C.7.  The  mean  of  the  sixth  and  seventh  summations  is  zero  by  Eq. 
C.8.  Most  of  the  terms  in  the  expansion  of  the  third  summation  are  zero  by  Eq.  C.4,  and  most 
of  the  terms  of  the  eighth  by  Eq.  C.9.  The  non-zero  portion  of  the  third  summation  is 
determined  by  the  application  of  Eq.  C.5,  and  the  non-zero  portion  of  the  eighth  summation 
by  Eq.  C.10.  The  application  of  Eqs.  C.S  and  C.10  to  these  terms  has  the  effect  of  converting 
the  third  summation  of  Eq.  28  from  a  double  to  a  single  sum,  converting  the  eighth 
summation  from  four-fold  to  two-fold,  and  eliminating  all  the  exponential  factors.  Combining 
some  of  the  terms  remaining  algebraically,  and  using  Eqs.  25  and  26  as  well  as  the  result  that 
the  wave-wave  interaction  function  R^Otj.Kj)  is  zero  (appendix  A),  results  in  the  following 
Fourier  sum  approximation  to  the  statistical  mean  value  of  Eq.  28: 

(?(t))  =  ^  1  ftt(<cj,gj.*j,Q'j)  Tt(«j)  Fw(<Cj)  A?Cj 

i 

+  2  +  +  +  +  Q^(*j,Km)  F„ («j)  FK(*m)  A*, A*ro 

J  **» 

^  1  —  —  ^T^J^m)  f  !»(*“))  A«j  A«m 

J  *»» 

(29) 
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The  compress  notation  in  this  result  the  wave  number  window  functions  have  been 
combined  with  the  wave-wave  interaction  functions  in  accordance  with  Eqs.  8  through  10,  and 
component  wave  number  intervals  are  represented  by  vector  wave  number  intervals.  Equation 
25  and  previous  definitions  of  the  wave-wave  response  functions  make  the  expression  purely 
real. 

Because  all  the  complex  random  amplitudes  have  been  eliminated,  the  range  of  the 
summations  may  be  allowed  to  increase  and  the  wave  number  intervals  allowed  to  decrease  so 
that  the  limiting  form  of  Eq.  29  involves  ordinary  integrals  as  follows: 

(t(t))  =  f...fni(Kl,(ri,Ki,al)T%(Kl)rn(Kl)dKl 

+ f...fll((Ki+Ke.tTt+<Ta.K1+KB,0t+ae)  Q5(*i,*8)  Fm(«i)  F»(*b)  d5gl  d*B 

+ y\../Rt(*1-*8,o1-a8,«1-*B.o1-ffa)  Rt(«i.*s)  Fm(*i)  Fh(*8)  d*i  d*a  (30) 

Apart  from  the  wave  number  window  Tj(2)  the  first  term  of  Eq.  30  is  exactly  the  same  as 
the  estimator  for  the  mean  developed  earlier  [3]  for  the  linear  wave  field  case.  Since  the  wave 
number  window  has  been  assumed  to  be  unity  up  to  some  large  wave  number,  its  inclusion  in 
the  first  term  makes  no  practical  change  in  the  estimator  because  the  nominal  wave  number 
variance  spectrum  FN(K)  generally  tends  toward  zero  for  large  wave  number. 

Comparison  with  Eq.  20  and  the  discussion  of  Eq.  28  show  that  the  last  two  terms  of  Eq. 
30  are  the  mean  quadratic  ship  system  response  to  the  true  wave  field  components  produced  by 
wave-wave  interaction.  The  form  is  very  like  that  of  the  true  variance,  Eq.  18.  Were  it  not  for 
the  sum  and  difference  frequency  arguments  in  the  ship  system  response  function,  it  would  be 
possible  to  change  variables  in  Eq.  30  so  as  to  compact  the  expression. 

On  the  assumption  that  the  order  of  magnitude  of  the  terms  in  Eq.  30  is  determined  by 
the  magnitude  of  the  wave  excitation,  the  first  term  is  second  order  and  the  last  two  terms  are 
fourth  order.  In  effect,  if  the  representation,  Eq.  28,  is  assumed  to  be  correct  through  fourth 
order,  the  modification  of  the  nominal  wave  field  by  wave-wave  interactions  influences 
estimates  of  the  mean  response  only  by  the  addition  of  terms  which  are  two  orders  of 
magnitude  smaller  than  the  leading  term. 

STATISTICS  OF  RESPONSE  TO  SECOND  ORDER 

In  the  discussion  of  the  spectral  representation,  Eq.  28,  it  was  noted  that  this 
representation  has  been  developed  on  the  assumption  that  the  truncated  functional  expansions 
of  both  the  true  wave  field  and  the  response  are  exact.  Had  the  functional  expansions  been 
truncated  at  third  degree  nonlinearities,  for  example,  two  additional  third  degree  contributions 
would  appear  in  Eq.  28.  These  contributions  would  correspond  to  the  linear  response  to  third 
order  wave-wave  interactions  and  a  response  of  third  degree  to  the  nominal  linear  wave  field. 
There  would  also  appear  contributions  to  the  terms  of  fourth  degree,  as  well  as  new  fifth  and 
sixth  degree  terms.  If  the  order  of  magnitude  of  the  various  contributions  is  assumed  to  be 
based  on  the  wave  components  as  outlined  previously,  Eq.  28  includes  third  and  fourth  order 
contributions  which  are  quite  possibly  incomplete,  and  thus  inconsistent. 
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Inconsistency  through  nonlinearities  of  second  degree  (order)  [3]  led  to  the  present  work. 

By  exactly  the  same  arguments  the  present  results  through  Eq.  30  are  inconsistent  with  respect 
to  third  and  higher  order  responses. 

Thus,  to  develop  response  statistics  which  are  consistent  to  second  order  it  is  necessary  to 
disregard  the  terms  of  Eq.  28  of  third*  and  higher  degrees.  Accepting  this  premise,  and  defining 
|fc)  as  the  Fourier  sum  approximation  to  the  spectral  representation  of  the  response  through 
second  order  results  in 

?(t)  -  y)  Bxp[jgjt]  ?(gjta j)  | 

J 

+  s’S  £  J)  W*(m)  Bxp[i(o j+om)t]  Q(«j,Oj,icm,am)  \ 

+  g£  £ W(m) ■»p[i(rr,-om)t)  R(7cJ,crJ,icm.am)  j  (31) 

J  m  . 

where  the  products  of  wave  number  windows  and  response  functions  have  been  consolidated  as 
follows: 


T(2|.Oj)  -  TH(#cJ)Tt(*J,aJ)  (32) 

Q(<j*<rJ,«m,o,m)  *  T„(XJ)TI((*m)x 

[Q{(cJtXm)  +Qt(xjfff  JfEm,am)]  (33) 

R(*J.ffJ.*m.am)  -  Tlf(«J)TII(«m)x 

s*m)  Tj(rcj  Km,(Tj  u jj,)  +R|(iCj((7 jirmiOn)) J  (34) 

It  will  be  noted  that  the  functions  defined  in  Eqs.  32  through  34  are  effective  linear  and 
quadratic  ship  system  response  functions  relating  the  response  to  the  nominal  linear  wave  field 
components.  When  the  wave  number  window  is  assumed  real  and  equal  to  unity  over  the  range 
of  wave  number  of  interest,  the  effective  linear  function,  Eq.  32,  is  the  same  as  the  linear 
response  function  defined  in  Eq.  22.  Apart  from  notation  it  is  the  same  as  the  linear  response 
function  previously  developed  [3].  The  quadratic  response  functions,  Eqs.  33  and  34,  have  the 
same  symmetries  as  the  response  functions,  Eq.  25,  by  virtue  of  the  symmetry  of  the  wave- 
wave  interaction  functions  and  the  definition,  Eq.  22.  In  addition,  because  of  the  zero  mean 
property  of  the  wave-wave  difference  frequency  response  function, 


RUj.Oj.cj.Oj)  =  T^(«j)  Rf(jCj,Oj,#cj,<7j) 


(35) 


In  effect,  the  function  determining  the  mean  value  of  the  effective  response  to  a  true  wave 
component  is  the  same  as  the  corresponding  function  for  the  response  to  a  linear  wave 
component. 
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As  would  be  expected,  the  spectral  representation,  Eq.  31,  is  analogous  to  the  second 
order  system  expansions  of  Papanikolaou  [14],  in  which  the  contributions  of  second  order 
nonlinearities  include  both  pure  second  order  hydrodynamic  response  and  interactions  of 
second  order.  In  the  present  case  the  grouping  of  the  origins  of  the  response  is  less  specific. 

(No  attempt  has  been  made  here  to  analyze  the  hydrodynamic  response  problem  other  than  to 
note  that  the  “true  wave  field”  could  just  as  well  have  been  represented  by  “true  wave 
potential”.) 

Equations  33  and  34  explain  what  Hasselmann  [1]  meant  by  the  assertion  that  the 
corrections  to  the  response  due  to  wave-wave  interactions  could  be  absorbed  into  the  quadratic 
response  functions.  In  a  consistent  expansion  correct  to  second  order  in  the  true  wave  field  the 
wave-wave  interaction  functions  can  be  combined  with  the  ship  system  response  functions  with 
advantage.  In  a  higher  order  expansion  the  second  order  wave-wave  interactions  also  contribute 
to  the  higher  order  response  functions,  so  that  the  details  of  the  “absorbtion”  are  relative  to 
the  order  of  the  expansion. 

As  just  noted,  the  symmetries  of  the  effective  functions,  Eqs.  32  through  34,  are  the  same 
as  the  ship  system  response  functions,  Eqs.  22  through  24.  Both  these  sets  of  functions  are  the 
same  as  the  linear  and  quadratic  response  functions  developed  earlier  [3]  in  the  sense  that  they 
represent  the  normalized  complex  response  of  the  system  to  linear  components  of  the  wave 
field.  The  spectral  representation  through  second  order  of  the  response  to  the  true  wave  field, 
Eq.  31,  is  identical  in  form  to  the  earlier  spectral  representation  of  the  response  to  the  linear 
field  (Eq.  34  of  reference  [3]).  The  difference  is  that  the  linear  and  quadratic  response 
functions  of  Eq.  31  represent  the  response  to  both  the  linear  field  and  the  second  order 
components  produced  by  wave-wave  interactions.  The  subsequent  derivations  of  statistics  of 
response  in  the  earlier  work  [3]  involved  only  the  symmetries  of  the  response  functions.  Thus  it 
is  pointless  to  continue  here  with  derivations  of  the  various  statistics.  The  results  will  be  the 
same  as  those  of  the  earlier  work  [3]  when  the  effective  response  functions,  Eqs.  32  through  34 
are  substituted  for  the  corresponding  response  functions  of  reference  [3]. 

With  respect  to  the  most  practically  important  statistic,  the  expected  value  of  response,  the 
above  result,  with  Eq.  35,  means  that  to  second  order  the  estimator  for  the  mean  response  to 
the  true  wave  field  is  exactly  the  same  as  the  estimator  for  the  mean  response  to  the  linear 
field.  To  this  order  the  mean  value  of  response  is  not  affected  by  wave-wave  interactions. 

To  the  extent  that  the  effective  response  functions,  Eqs.  32  through  34,  are  obtainable,  the 
other  statistics  derived  earlier  [3]  are  also.  Because  the  wave  number  windows  in  Eqs.  32 
through  34  have  consistently  been  assumed  to  be  unity  in  the  wave  number  range  of 
importance,  they  constitute  no  particular  problem  —  in  fact,  they  formally  represent  only  the 
finite  truncation  of  wave  number  interval  which  would  be  necessary  in  a  practical  calculation. 
On  this  basis  Eq.  32  and  the  second  terms  of  Eqs.  33  and  34  also  constitute  no  new  problems. 

The  wave-wave  interaction  functions  Q{(xj,xm)  and  R^(xj,xm)  appearing  in  Eqs.  33  and  34 
have  been  defined  explicitly  for  deep  water.  The  linear  response  functions  of  the  form 
T4(xj±xm,Oj±orn)  which  appear  in  Eqs.  33  and  34  present  a  slightly  different  problem  from 
the  ship  motions  perspective.  The  argument  Xj±3cm  defines  a  wave  which  is  in  general  of 
different  length  and  direction  than  either  of  the  superimposed  linear  waves.  However,  the 
argument  oj±om  denotes  that  the  response  is  at  a  frequency  which  is  the  sum  or  difference  of 
the  encounter  frequencies  of  the  two  linear  waves.  The  phase  speed  of  the  second  order 
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interaction  component  is  not  determined  by  the  free  wave  dispersion  relationship.  This  result  is 
consistent  with  the  result  that  the  second  order  component  of  a  Stokes  wave  has  the  same 
phase  speed  as  the  first  order  component,  although  the  wave  number  is  twice  that  of  the  first 
order  component.  Obtaining  the  required  linear  response  function  of  the  ship  motion  system,  if 
it  is  needed,  will  thus  require  the  treatment  of  a  somewhat  different  potential  flow  situation 
than  is  ordinary.  The  level  of  effort,  however,  is  likely  to  be  much  smaller  than  the  effort 
required  for  the  quadratic  responses  to  the  linear  waves. 

CONCLUSIONS 

The  present  objectives  included  a)  the  development  of  a  model  for  the  statistics  of 
quadratic  ship  system  response  in  short-crested  seas  which  is  consistent  in  that  quadratic 
nonlinearities  in  the  wave  field  are  taken  into  account,  b)  a  clarification  of  the  assumptions 
being  made,  c)  an  assessment  of  what  is  involved  in  carrying  out  computations  which  include 
the  wave  field  nonlinearities,  and  d)  an  assessment  of  the  influence  of  the  wave  field 
nonlinearities  upon  the  most  important  statistic  associated  with  quadratic  response,  the  mean 
value. 

In  this  development  the  lead  of  Hasselmann  was  followed  and  amplified.  If  the  statistical 
representation  of  the  response  to  the  linear  part  as  well  as  to  the  quadratic  nonlinearities  of  a 
deep-water  wave  field  is  carried  out  strictly  through  contributions  which  may  be  considered 
second  order  in  wave  excitation,  and  if  higher  order  response  is  neglected,  the  statistical 
estimating  formulae  for  the  quadratic  response  to  short-crested  seas  in  which  quadratic 
nonlinearities  are  included  are  identical  to  those  previously  developed  for  the  case  in  which  no 
wave  nonlinearities  are  taken  into  account. 

The  only  difference  is  that  the  quadratic  sum  and  difference  frequency  response  functions 
required  by  the  previous  work  are  modified  by  a  term  which  amounts  to  the  linear  ship  system 
response  to  those  quadratic  wave  field  components  produced  by  the  wave  field  nonlinearities. 
The  main  assumption  producing  this  result  is  that  the  magnitude  of  the  nonlinear  ship  system 
response  follows  the  magnitude  of  the  wave  excitation,  and  thus  that  nonlinear  ship  system 
response  to  quadratic  wave  field  components  is  of  higher  order  and  may  reasonably  be 
neglected.  The  alternative  to  this  assumption  involves  the  development  of  a  consistent  higher 
order  model  which  has  not  been  attempted  here. 

In  so  far  as  practical  computation  of  response  functions  is  concerned,  the  conceptual  idea 
is  the  same  as  in  previous  work.  The  linear  response  function  is  exactly  the  same.  The  present 
development  requires  that  the  response  through  second  order  terms  in  excitation  be  computed 
for  pairs  of  waves,  just  as  in  the  previous  work.  The  difference  in  the  context  of  potential  flow 
hydrodynamics  is  that  the  assumed  wave  potential  for  the  wave  pair  must  include  the  potential 
due  to  the  quadratic  wave-wave  interactions.  The  net  effect,  after  the  computed  responses  in  sum 
and  difference  encounter  frequencies  have  been  identified,  will  be  to  add  the  linear  response  to 
the  wave  field  nonlinearities  noted  above.  This  additional  contribution  is  a  slightly  different 
potential  flow  problem  because  the  wave-wave  interaction  components  are  not  free  waves.  On 
the  whole,  however,  if  the  addition  of  this  effect  is  required,  it  could  easily  involve  less  work 
than  is  involved  in  the  computation  of  quadratic  response  to  the  linear  components  of  the 
wave  field. 


The  most  important  statistic  in  practice  involving  quadratic  response  is  the  expectation,  or 
mean  value.  In  this  special  case  the  extension  just  outlined  is  not  necessary.  To  second  order, 
the  mean  response  of  a  linear  plus  quadratic  ship  system  to  a  short-crested  wave  Held  is  not 
affected  by  wave  field  nonlinearities.  The  previously  developed  approach  appears  to  apply 
equally  well  when  quadratic  wave-wave  interactions  are  taken  into  account.  The  present 
development  suggests  that  the  mean  value  estimator  may  be  robust  through  third  order, 
although  a  consistent  development  through  this  order  must  be  made  for  confirmation. 
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APPENDIX  A 

QUADRATIC  WAVE-WAVE  INTERACTION  THEORY 

The  theory  of  quadratic  wave-wave  interactions  has  been  developed  in  a  number  of  ways. 
The  treatment  here  summarizes  that  of  Longuet-Higgins  [9],  and  will  be  specialized  to  the  deep 
water  case.  To  first  order,  the  potential  solutions  for  small  amplitude  gravity  waves 
superimpose.  The  waves  propagate  independently  and  without  interaction.  To  second  order  the 
waves  interact,  and  the  interaction  produces  a  small,  bounded  modification  to  the  wave 
motion.  It  is  the  second  order  modifications  which  are  of  interest  here,  and  to  study  these  it  is 
sufficient  to  consider  the  interactions  between  pairs  of  waves. 

To  fix  notation  the  parameters  describing  a  single  small-amplitude  gravity  wave  will  be 
taken  as  follows: 


X  = 

wave  length 

Cl.  = 

wave  frequency 

K  = 

vector  wave  number 

with  components; 

*X  = 

*cos(/i) 

*Y  = 

#c  sin(ji) 

M  = 

wave  direction 

= 

tan_,(«Y//cx) 

*  = 

scalar  wave  number 

= 

2tr/X  =  |*| 

= 

«a/g  for  deep  water 

where  the  deep-water  dispersion  relationship  just  expressed  will  be  assumed  hereafter,  and  a 
fixed  X-Y  coordinate  system  in  the  undisturbed  water  surface  with  the  Z  coordinate  positive 
upward  is  implied. 

The  derivation  of  Longuet-Higgins  involves  the  postulation  of  a  velocity  potential 
satisfying  the  Laplace  equation  and  the  free  surface  boundary  conditions.  The  result  is  three 
equations  in  the  unknown  potential,  the  velocity,  and  the  free  surface  elevation.  These 
equations  are  expanded  in  a  Taylor  series  about  the  mean  water  position,  and  formal 
perturbation  series  are  assumed  for  the  potential,  velocity,  and  surface  elevation  of  two 
intersecting  trains  of  regular  waves.  Substitution  of  the  series,  gathering  terms  of  first  and 
second  order,  with  the  requirement  that  the  potential  satisfy  the  Laplace  equation  results  in  a 
solution  for  the  total  velocity  potential  of  the  two  intersecting  waves  (correct  to  second  order 
for  deep  water)  in  the  form 

^  —  ^01  +  ^  10  ^08  ^ 20  +  (A.l) 

In  this  result  the  first  and  second  terms  are  the  well-known  potentials  for  two  first-order 
regular  waves.  The  third  and  fourth  terms  represent  quadratic  interactions  of  each  first-order 
wave  with  itself.  Both  are  functions  of  time  only,  corresponding  to  the  vanishing  of  second- 
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order  velocities  in  a  single  first-order  wave.  The  last  term  represents  the  potential  resulting 
from  second-order  interactions. 

To  be  specific,  the  two  intersecting  first-order  deep-water  wave  trains  may  be  assumed  in 
the  form 

Cio  =  *iCoa(*1  •  X  — Wjt  +  Ci) 

? 01  =  »8cos(*8  •  X  -  w8t  +  e8)  (A.2) 

where  X  represents  the  vector  position  in  fixed  coordinates  and  is  an  arbitrary  phase  angle. 
The  corresponding  first-order  potentials  may  be  taken  in  the  form 

io  =  ai«i<cr1Exp[*iZ]  ainf*!  •  X  -  Wit -J-Ci) 

V»oi  =  a8“a*81Exp[*8Z]  sin(«8*  X- wat +  ea)  (A.3) 

With  these  representations  and  the  result  for  the  interaction  potential  from  Longuet- 
Higgins  [9]  the  total  velocity  potential  for  deep  water  of  two  intersecting  wave  trains  may  be 
represented  by 

<p  =  a1«1/Ci1Exp[#c1Z]  ainf*!  •  X  -  Wjt  +  tj) 

+  a8wax;a1E*p[x;2Z]  sin(ic8  •  X  -  w8t  +  e8) 

+  at  a8  A(7Ci  ,*8)  Exp[|«i-«8|Z]  x 

sinj(xt  -  *8)  •  X  -  (wi  -  wa)t  +  (et  -  e8) \ 

-  aj  a8  B(xi  ,kb)  Exp[|x,4xg|  Z]  x 

sin|(«i  +  Ka)  •  X  -  («i  +  wa)t  +  (en  +  e2)|  (A.4) 

The  third  and  fourth  terms  of  Eq.  A.4  are  the  second-order  wave-wave  interaction  potential, 
and  the  as  yet  undefined  functions  in  these  terms  may  be  written 

SoiiWgfwi-wg)  coss[(/i .-a*2)/2] 

A(*i  ,kb)  =  - - - - - - - 

(Wi-W8)8  -  gl/Ci-Kgl  (A. 5) 

EWiWefwi+wg)  ain8[(A*1-A*B)/2] 

B(/c1(/c2)  =  - - - 

(wi  +  wg)  -g|«i  +  #68|  (A.6) 

With  the  velocity  potential  defined  the  total  surface  elevations  of  the  two  intersecting  wave 
trains  may  be  derived  from  the  Longuet-Higgins  analysis  [9]  and  are  written  as  follows: 
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*  •1co«(21  •X-a>1t  +  c1) 


+  i|Coi(Z|*X-»|t-t-C|) 
e»* 

+  »*  — coi|8(*1‘f)-8»1t  +  8cl| 
o»*  _ 

+  ag  “  coa|8(£a*  X)  -  2wgt  +  8cgj 

+  *1 »»  9*(*i  .*•)  ®oi j(*i  +  *g)  •  X  -  (Wj  +  «g)t  +  («!  +  eg)j 
+  *1  *a  Rt(*i  •*»)  co8|(*i  -  Eg)  •  X  -  (wj  -  «g)t  +  (ex  -  eg)|  (A.7) 

As  in  the  case  of  the  potential,  two  new  functions  appear  and  are  defined  as  follows: 


QfUi  .*«) 


Rf(*l  .*») 


7  +  l  +  > 

/  7*~ainBr 

i<n*rC 

7  +  l-> 

J  7a-ain,r 

■  1(1  I  / 

(A.8) 

7-l+> 

J-r*~co»Mr 

7 -l-\ 

]y*~co»*r 

vUb  a  r 

(A.9) 

where:  <un  =  ^kjg 

=  L( 

2  \  «B  «1  / 

r  -  (/*,-Mg)/2 

(A.  10) 

Special  limiting  cases  of  Eqs.  A.8  and  A.9  are  as  follows: 

=  wf/g 

=  0 

(A.l  1) 

The  result  for  the  corrected  wave  field,  Eq.  A.7,  is  very  close  in  form  to  that  obtained  [3] 
for  the  response  of  a  linear  plus  quadratic  system  to  two  superimposed  sinusoids.  If  the 
position  in  the  wave  field  is  taken  to  be  fixed,  the  expression  involves  six  frequencies:  those  of 
the  two  small  amplitude  waves,  their  second  harmonics,  and  the  sum  and  difference 
frequencies.  Missing  is  the  shift  in  the  mean  obtained  for  the  linear  plus  quadratic  system. 
However,  terms  corresponding  to  a  shift  in  the  mean  would  involve  the  special  case  of 

indicated  in  the  second  of  Eqs.  A.  11,  so  that  Eq.  A.7  could  be  put  in  the  form  of 
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the  response  of  •  liner  phis  quadratic  system.  Essentially,  the  analysis  says  that  there  is  no 
shift  hi  the  mean  of  the  corrected  wave  field  (as  might  be  expected  from  continuity 
considerations  for  incompressible  flow).  The  functions,  Eqs.  A. 8  and  A.9,  are  real  and 
symmetric  in  their  arguments,  and  thus  they  satisfy  the  requirements  for  quadratic  frequency 
response  functions  noted  in  Dalzdl  [3],  These  considerations  along  with  the  position  of  the 
functions  in  Eq.  A.7,  suggest  that  they  may  be  regarded  as  quadratic  frequency  response 
functions  which  relate  the  sum  and  difference  frequency  components  of  the  true  wave  field  to 
the  interactions  of  the  two  superimposed  small  amplitude  waves.  The  first  of  Eqs.  A.l  1  serves 
to  put  the  coefficients  of  the  second  harmonics  of  Eq.  A.7  into  the  form  of  the  response  of 
the  linear  plus  quadratic  system.  Eq.  A.7  contains  no  explicit  notation  for  the  linear  frequency 
response  function,  which  thus  may  be  regarded  as  unity. 

The  preceding  arguments  indicate  that  the  true  wave  Field  (to  second  order)  may  be 
regarded  as  the  response  of  a  linear  phis  quadratic  system  to  a  Active  “linear”  small  amplitude 
wave  field.  The  quadratic  system  is  defined  by  a  linear  frequency  response  function  of  unity 
and  quadratic  frequency  response  functions  defined  by  Eqs.  A.8  and  A.9.  Although  it  has  not 
been  done  explicitly  here,  the  wave  potential,  Eq.  A.4,  may  also  be  interpreted  in  this  way. 
These  observations  form  the  hydrodynamic  background  for  Hassdmann’s  approach  [1]  to  the 
consistency  problem  mentioned  in  the  main  body  of  the  report. 

Had  the  Longuet-Higgins  derivation  been  for  the  true  wave  potential  and  elevations  (to 
second  order)  corresponding  to  the  superposition  of  N  independent  small  amplitude  waves,  the 
perturbation  expansion  would  have  involved  the  interactions  of  all  possible  pairs  of  waves. 
These  interactions  would  have  given  rise  to  a  total  potential  of  the  form 


(A.  12) 


The  summations  on  pairs  of  waves  carries  through  to  the  elevations.  With  a  manipulation 
of  Eq.  A.7,  and  a  use  of  Eqs.  A.ll,  the  true  wave  field  corresponding  to  the  summation  of  N 
small  amplitude  waves  may  be  put  in  the  form 


N 

C(x,t)  =  yt  J?g|bjExp[iCj  •  I  -  ia>jt  ]| 

J 

a  a 

Exp[i(«j  +  **)•  X  —  i(ci)j  +  «k)t]j 

*  N  * 

J  a  _ 

Exp[i(«j-*k)*  X  -  i(&>j  -  «n)t]j  (A.  13) 

where  bj  =  ajExp[Uj] 


With  the  position  vector  fixed,  this  result  is  nearly  identical  to  that  for  the  response  of  a 
linear  plus  quadratic  system  to  the  sum  of  N  small  amplitude  waves. 
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APPENDIX  1 

FUNCTIONAL  EEFBESENTATION  OF  THE  WAVE  FIELD 

The  hydrodynamic  arguments  suggest  that  the  true  wave  field  may  be  taken  to  be  a 
functional  of  the  “linear”,  small  amplitude  Add,  at  least  to  second  order.  In  accordance  with 
the  previous  work  involving  system  response  [3],  if  the  linear  field  is  continuous  and 
differentiable  and  vanishes  outside  some  finite  region  in  the  space  of  all  admissible  linear 
fields,  or  if  the  field  is  square  integrable,  the  true  wave  field  may  be  represented  as  a 
functional  Taylor  expansion  in  the  neighborhood  of  zero  wave  field.  With  a  further 
assumption  to  be  noted,  the  Taylor  expansion  through  the  second  order  term  may  be  written  as 
follows: 

fT(t.XvY)  -  fo 

+///* *xt  1)  * .X“*i  .T-y jJdtj  dxt  dy  x 

+/"  /*s<t i.*i.yi.*s.xa.yg)  * 

fL(t-tl,X-*1,Y-y1)fL(t-t,,X-x,.T-yg)  x 

dtj  dxt  dy  t  dts  dx,dyt  (B.l) 

where  the  integrals  are  infinite.  The  kernels'  within  the  integrals  in  Eq.  B.l  are  functions  of 
time  and  space  differences  and  contain  the  physical  relationships  between  the  linear  and  true 
wave  fields.  The  variables  t,  X,  and  Y  denote  absolute  time  and  position  in  a  fixed  frame  of 
reference.  As  in  the  previous  work  it  has  been  assumed  that  the  integrands  differ  from  zero 
only  in  a  finite  range  of  the  dummy  arguments. 

The  difference  between  Eq._  B.'l  and  the  corresponding  representation  [3]  of  the  response 
should  be  noted.  The  response  of  the  ship  system  is  supposed  invariant  with  absolute  position, 
and  in  the  representation  of  the  response  this  restriction  was  taken  into  account  in  the 
formulation  for  the  encounter  domain.  In  effect,  had  the  response  been  presumed  a  function 
of  absolute  position,  the  form  of  the  final  functional  representation  would  have  been  of 
exactly  the  same  form  as  Eq.  B.  1 

The  hydrodynamics  of  the  problem  noted  in  Appendix  A  allows  an  immediate 
simplification  of  the  representation.  First,  note  that  in  Eq.  B.l,  when  the  linear  field 
approaches  zero,  there  is  left  the  zero’th  order  term  go,  and  this  term  may  be  taken  equal  to 
zero  on  the  basis  that  the  true  wave  field  must  approach  the  linear  field  as  the  latter 
approaches  zero.  The  hydrodynamic  analysis  also  shows  that  there  is  no  linear  transformation 
of  the  linear  field.  Thus  the  kernel  g](t|,X|,y|)  of  the  second  term  of  Eq.  B.l  may  be  regarded 
as  a  defta  function  in  the  dummy  variables.  The  net  result  of  these  considerations  is  that  the 
*  functional  representation  of  the  true  wave  field  becomes 


fT(t,X.Y)  -  CL(t,X.Y) 

CL(t-t,  .X-Xi  .Y-,1)CL(t-ti.X-x..Y-yi)  x 

dt,  d*t  dy,  dtadxadyg  (B.2) 

As  in  the  development  of  the  interpretation  of  the  lime  domain  kernels  of  response,  the 
quadratic  kernel  of  Eq.  B.2  may  be  interpreted  by  assuming  that  the  linear  wave  field  is  the 
superposition  of  two  small  amplitude  regular  waves,  which  is  the  same  as  Eq.  1  with  N  *  2. 

The  steady  state  true  wave  field  is  obtained  by  substitution  into  Eq.  B.2.  The  result  of  this 
operation  is  the  same  as  Eq.  2,  or  Eq.  A.  13,  with  N  =  2  if  the  following  transforms  hold: 

Qf(*i.*«)  =  /••/*,( t1.xt.j1.t,.xa.ya)  x 

Exp[-i(czaxa  +  cTaya-a>ata)]  dtj.-.dy,  (B.3) 

=  /  ./«B(t i,xl.yl.ta,xa.ya)  x 

E*P[-i(*xi*i  +  *T|y1-«»it|)]  x 

Exp[i(*zaxa+#cTaya-wata)]  dt,...dyf  (B.4) 

In  these  transforms  it  is  implicitly  assumed  that  frequency  and  wave  number  are  connected 
by  a  dispersion  relation.  Thus  the  sum  and  difference  response  functions  are  denoted  only  as 
functions  of  vector  wave  numbers,  although  this  is  not  mathematically  essential. 

To  ensure  that  Eqs.  B.3  and  B.4  are  valid  representations  the  integral  of  the  absolute 
value  of  the  time-space  kernel  in  the  transforms  must  exist,  and  the  kernel  must  be  continuous 
[11].  Additionally,  if  the  inverse  transforms  are  valid,  the  kernel  must  be  square  intcgrable, 
and  the  same  continuity  and  integrability  conditions  must  hold  for  the  inverse  functions  of 
wave  number. 

These  last  considerations  pose  a  serious  practical  problem  in  view  of  the  results  of  the 
hydrodynamic  theory.  Appendix  A.  The  hydrodynamic  functions  corresponding  to  the  left- 
hand  sides  of  Eqs.  B.3  and  B.4  are  defined  in  Eqs.  A. 8  through  A.ll.  As  the  magnitude  of  the 
vector  wave  number  arguments  tends  to  infinity,  the  value  of  the  functions  does  also.  Thus  the 
functions  defined  by  Eqs.  A.8  and  A.9  do  not  satisfy  the  conditions  required  for  the  validity 
of  Eqs.  B.3  and  B.4.  As  a  practical  example,  were  it  required  to  generate  the  time-space  kernel 
by  transforming  the  response  functions,  there  would  be  an  immediate  problem  because  the 
limits  of  the  integrals  of  Eqs.  B.3  and  B.4  as  well  as  those  of  their  inverses  are  infinite. 

Hineno  [7]  ran  into  this  problem  in  conjunction  with  a  simulation  of  long-crested  waves  and 
solved  it  numerically  with  a  truncation  technique  which  is  mathematically  similar  to  that  which 
will  be  developed  next. 
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One  way  out  of  such  a  problem  is  analogous  to  the  “white  aoise”  excitation  approach 
often  used  to  mitigate  similar  problems  in  the  time-only  case.  Essentially,  in  the  time-only  case 
white  noise  is  operated  upon  with  a  linear  transformation  to  produce  an  excitation  to  the 
system  of  interest,  and  the  response  functions  connecting  the  noise  to  the  system  response  are 
considered  in  bcu  of  those  connecting  the  actual  excitation  to  the  response.  In  the  time-only 
case  the  implementation  is  through  well-known  relationships  developed  for  the  results  of 
cascading  Voherra  series. 

Thus  to  pursue  this  direction  in  the  present  case  it  is  necessary  to  hypothesize  a  noise  or 
nominal  field  £^,5?)  in  which  time  and  position  vectors  have  the  same  meaning  as  before.  The 
cascade  of  fields  may  be  diagrammed  as  follows: 


*  fl(.) 

gt(.) 

«.<)_ 

The  nominal  field  is  presumed  to  be  operated  upon  by  a  linear  system  ft(.)  in  the  time/space 
domain  to  produce  the  linear  small  amplitude  wave  field  £L(t,X)  which  is  defined  as  before. 
The  linear  small  amplitude  field  is  acted  upon  by  a  linear  plus  quadratic  system  with  kernels 
g,(.)  and  g2<.)  to  result  in  the  true  wave  field  £-j-(t,X!).  The  last  functional  relation  is  the  same 
as  Eq.  B.  I  with  the  zero’th  term  omitted. 

The  relationship  desired  between  the  nominal  field  and  the  true  wave  field  may  be 
diagrammed  as  follows: 


h,(.) 

CTU.X) 

h.(.) 

where  the  kernels  h,(.)  and  hjf.)  incorporate  the  relationships  between  the  nominal  and  linear 
fields  and  those  between  the  linear  and  true  wave  fields. 

If  the  zcro'th  order  term  is  omitted,  the  corresponding  time/space  domain  Taylor 
expansion  of  the  true  wave  field  in  terms  of  the  nominal  field  has  the  same  form  as  Eq.  B.l, 

CT(  t.  X.  Y )  -  /// hl(t1.xl.yl)CN(t-tl.X-xl.r-y|)dtldx1dy| 

♦/  /M11  .*1  .y,  .ts.*«.y,)  x 

C„(t  t,,X  «„T  y,K1,(t-t,,X-xg.Y-y#)  x 

dt,  dx,  dy,  dt,dx,dyB  (B.5) 

On  the  assumption  that  the  kernels  in  Eq.  B.5  are  sufficiently  well  behaved,  the  corresponding 
linear  and  quadratic  response  functions  relating  nominal  and  true  wave  fields  may  be  written 
as  follows: 
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Tt(1)  -  |.Z|.yt)  Bxp[-i(*xx,+*TyI-»t1)]  dt1dx1dy1 


(B.6) 


-  /•  / h.(t1,x1,y1,ti.x„yt)  x 

B*Pl-i(*xi*i+*T1y1-«|t1)]  x 

Exp[-i(KXBxg+cTays-wgtg)]  dtt...dy,  (B.7) 

Rr(*i.i«)  =  /  /hg(t1,x1.y1.tg.xg.y§)  x 

K*p[-i(*I1*i+*Tiy1'wit1)]  x 

Exp[i(*laxa+*Tiyg-ft>gtg)]  dti ...dyg  (B.8) 

To  complete  the  definition!  required  for  the  cascade  problem  the  assumed  linear 
transformation  between  noise  and  linear  wave  fields  must  be  defined  and  may  be  written  in  the 
time/space  domain  as  follows: 

tL(t.X.Y)  =  ///f1(t1,x1.y1)fI((t-t1,X-x1.Y-y1)dt1dxIdy1  (B.9) 

The  linear  response  function  corresponding  to  the  kernel  of  Eq.  B.9  will  be  written 

T*(«)  =  /••/ fiUt.Xs.yj)  Expl-KxxXj  +  Xyyj-wt,)]  dt1dx,dy1  (B  IO) 

Now  if  the  model,  Eq.  B.9,  for  the  linear  field  is  substituted  into  the  initial  expansion  of 
the  true  field,  Eq.  B.2,  the  dummy  time  and  space  variables  are  renumbered,  and  a  simple 
variable  transformation  is  made  (noting  that  all  integrals  are  infinite),  the  result  is  sufficiently 
similar  to  the  expansion,  Eq.  B.5,  that  the  relations  between  the  kernels  of  Eqs.  B.2,  B.9,  and 
B.5  are  developed.  These  turn  out  to  be 

hi(t,.xi.y,)  =  Mti'Xi.yj)  (B.ll) 

h«(li  .*1  .yt  .t*.x*.yg)  - 

/  /*g(ti-t3.xi-x3.y1-y3.ts-t«.x8-x«.yg-y4)  * 

f|(t3.*3.y3)  f  i(*«  .*«  .y  4)  dt3dx3...dy4  (B.12) 

The  relations  between  the  response  functions  of  the  cascade  and  those  relating  the  nominal 
and  true  wave  fields  are  developed  by  substituting  Eqs.  B.ll  and  B.12  into  the  definitions, 

Eqs.  B.6  through  B.9,  with  the  following  results: 

Tt(*)  =  T*(ic)  (B.13) 
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(B.  14) 


=  Qf(*1.*«)T,|(21)T.(ie,) 

Tjj(*j)  Tj|(X|)  (B.15) 

If  the  functions  defined  in  Eqs.  B<14  and  B.15  can  be  made  to  satisfy  the  square 
integrability  condition,  the  representation,  Eq.  B.5,  will  be  valid,  as  will  be  the  inverse 
transforms  of  Eqs.  B.7  and  B.8.  For  this  to  be  true,  the  kernel  of  Eq.  B.9  and  its  transform 
must  first  be  square  integrate.  Inspection  of  Eqs.  A.8  and  A.9  discloses  that  under  the  worst 
circumstance  (one  wave  number  very  much  greater  than  the  other)  the  values  of  the  wave-wave 
interaction  response  functions  are  nearly  proportional  to  the  greater  of  the  two  scalar  wave 
numbers.  Thus  if  TnOD  approaches  zero  at  a  rate  greater  than  1/|X|  for  large  wave  number, 
the  left-hand  sides  of  Eqs.  B.14  and  B.15  should  approach  zero  sufficiently  rapidly  to  satisfy 
the  integrability  condition. 

It  is  possible  to  interpret  the  linear  transformation  hypothesized  here  in  a  number  of  ways. 
If  the  vector  wave  number  spectrum  of  the  linear  field  may  be  put  in  the  form 
then  the  linear  transformation  could  be  called  a  “spectral  filter"  in  direct  analogy  to  the 
common  practice  in  the  time-only  case.  (Scalar  wave  number  spectra  of  ocean  waves  decay 
more  rapidly  with  wave  number  than  is  required  to  ensure  validity.)  Alternatively,  the 
transformation  might  be  viewed  as  an  analytical  “window"  which  limits  the  range  of  possible 
wave  interaction.  In  this  interpretation  TN(E)  could  be  chosen,  for  instance,  to  be  purely  real 
and  equal  to  unity  up  to  some  cutoff  wave  number  higher  than  the  range  of  interest,  after 
which  a  relatively  arbitrary  transition  to  zero  could  be  made.  The  present  development  is 
essentially  an  artifice  to  overcome  the  existence  problems  brought  about  by  the  hydrodynamic 
development  for  the  wave- wave  interactions. 

As  usual,  the  functions  defined  by  Eqs.  B.6  and  B.7  may  be  interpreted  by  assuming  that 
the  nominal  field  is  a  summation  of  cosinusoidal  waves  in  the  form 

N 

C„(t,X)  =  fle |cjExp[iiCj  •  X  -  i«jt  (B.  16) 

J 

where  Cj  =  ajExp[iCj] 

and  aj  and  cj  are  the  arbitrary  amplitudes  and  phases  of  the  N  cosinusoids.  When  Eq.  B.16  is 
substituted  into  the  model,  Eq.  B.5,  and  the  transform  definitions  are  taken  into  account, 
there  results  the  following  expression  for  the  true  wave  field: 
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X 

ff(t>X)  —  ^0{cj  Tn(Cj)  Exp[i*j  •  X  —  ici)jt]| 

J 

X  X 

+  z!2  2  *9ie* Ck  Qc  (*i  -*k)  t"(*i) 

J  *  . 

Exp[i(«j  +  *k)  •  X  -  i(«j  +  ct>k)t]| 

1  x  x 

+  2'2£/?eiC^Ck  Rf<*J’*k>Tn(*j)Tx  («k) 

J  k  _ 

Exp[i(«j  -  7ck)  •  X  -  i(Wj  -  a>k)t]j  (B.17) 

The  difference  between  the  results  of  the  hydrodynamic  analysis,  Appendix  A,  and  those 
just  completed  is  that  the  functions  defining  the  sum  and  difference  wave  interactions  are 
modified  by  the  linear  response  function  presumed  to  relate  the  nominal  and  linear  fields. 
Additionally,  the  linear  components  of  the  true  field  involve  an  arbitrary  linear  transformation 
not  present  in  the  development  of  Appendix  A. 
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APPENDIX  C 

SPECTRAL  REPRESENTATION  OF  THE  NOMINAL  WAVE  FIELD 

For  present  purposes  the  existence  of  a  nominal  underlying  linear  wave  field  {NfX.Y.t)  is 
assumed.  It  is  further  assumed  that  the  elevation  variance  of  the  nominal  wave  field  may  be 
represented  by 

•ft  ftft 

Var(f„)  =  f  f  P*(<cx  ,*T)d/cx  d/cT  (C.l) 

where  F^^xy)  is  a  vector  wave  number  spectrum. 

It  is  also  assumed  that  the  nominal  field  is  dispersive,  homogeneous  in  space,  stationary  in 
time,  and  Gaussian.  These  assumptions  together  imply  that  the  field  has  a  “spectral 
representation”.  Although  this  representation  is  often  written  in  the  form  of  a  Fourier-Stieltges 
integral,  it  is  convenient  to  work  with  the  corresponding  approximating  Fourier  sum. 
Accordingly,  the  approximating  Fourier  sum  of  the  spectral  representation  of  the  nominal  field 
may  be  written 

?„(X,Y,t)  =  ^/?e{w(j)Exp[i*XJX  +  i*TJY-ia>Jt]j  (C.2) 

J 

where  W(j)  is  a  complex  random  function  of  the  interval  AjcXj,AxYj.  This  function  has  the 
following  properties:  (<•>  denotes  the  ensemble  expected  value.) 


(W(j))  as  0 

( W(j)  W*  (k))  -  0  (for  J  *  k) 

{  W(J)  W*  ( j))  =  2F]l(icXj  ,kTj)  A#cZjA#cTj 
(w(j)W(k))  s  0 


(C.3) 

(C.4) 

(C.5) 

(C.6) 


Now  the  initial  assumption  on  the  nominal  wave  field  of  stationarity  in  time,  with  the 
Gaussian  process  assumption,  as  well  as  that  of  homogeneity  imply  that  the  random  functions 
W(j)  are  Gaussian  and  statistically  independent  for  different  wave  numbers  and  directions  of 
propagation.  Thus  from  the  properties  of  the  third  and  fourth  moments  of  Gaussian  variables 
the  following  additional  relations  arise: 

<W(j)W(k)W(m))  =  (w(j)W(k)W*(m))  s  0  (C.7) 

(w(J)W(k)W(m)W(n))  =  (w(j)W(k)W(m)W*(n))  =  0  (C.8) 

(w(j)W(k)W*(m)W*(n))  *  0  (C.9) 
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unless  (j.k  =  m,n  or  j,k  =  n,m)  in  which  case 
(w(j)W(k)W*(j)W*(k)) 

=  4Fn(/cXj  ,#cTj)FN(A;Xk  ,/cTk)  AkXjA«tj  A/cIkA#cYk 
=  8[FN(*Xj Yj)  A#eXjAcTj]8  (  if  k  =  j)  (C.10) 

The  complex  random  amplitudes  W(j)  are  assumed  to  be  Gausian  and  zero  mean.  Eqs. 

C.3  through  C.10  define  the  expectations  of  first  through  fourth  order  moments.  Under  this 
assumption  the  expectations  of  fifth  and  sixth  order  moments  are  also  known.  The  fifth  order 
expectations  are  identically  zero  because  of  the  Gaussian  assumption. 

The  sixth  order  expectation  of  Gaussian  variables  may  be  written  as  the  sum  of  fifteen 
products  of  three  .second  order  expectations.  The  pairing  of  variables  in  the  second  order 
expectation  products  corresponds  to  all  the  ways  in  which  the  six  variables  may  be  grouped 
into  three  pairs.  Noting  Eqs.  C.l  through  C.3,  the  only  way  the  second  order  expectation  can 
be  other  than  zero  is  when  one  of  the  terms  in  the  product  is  a  complex  conjugate.  Thus  the 
only  way  the  product  of  three  second  order  expectations  can  be  other  than  zero  is  that  each 
second  order  expectation  must  be  of  the  form  of  Eq.  C.5.  Accordingly,  exactly  three  of  the  six 
variables  in  the  sixth  order  expectation  must  be  complex  conjugates.  Otherwise  the  process  of 
grouping  the  six  variables  into  pairs  will  result  in  at  least  one  of  the  terms  in  each  of  the 
fifteen  triple  products  of  second  order  expectations  being  of  the  form  of  Eq.  C.6,  and  thus  the 
entire  expectation  will  be  zero. 

Writing  out  the  expansion  of  the  sixth  order  expectation  in  which  exactly  three  of  the 
variables  are  complex  conjugates,  and  applying  Eq.  C.2  results  in  the  following  definition  of 
the  only  ways  in  which  this  expectation  can  be  other  than  zero: 

(w(j)  W(k)  W(m)  W*(n)  W*(p)  W*(q)) 


)FN 

(« 

Xk  >'cTk)F 

'*Tm)  A/CXj 

A«Tj 

A*Tk 

If 

n 

*  j  • 

P  = 

k 

,  and 

q 

=  m 

» 

or 

n 

=  j  , 

P  = 

m 

,  and 

q 

=  k 

» 

or 

n 

=  k  , 

P  = 

m 

,  and 

q 

=  j 

> 

or 

n 

=  k  , 

P  = 

j 

,  and 

q 

=  m 

> 

or 

n 

=  m  , 

P  = 

j 

,  and 

q 

=  k 

» 

or 

n 

=  m  , 

P  = 

k 

,  and 

q 

=  J 

• 

(C.ll) 
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